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Xfy  f  ABSTRACT*  * 

U  f 

Experiences  learned  at  the  AFRL  on  injection  ofUcryogenic  jetymto  sub-  and  supercritical  conditions  using  shadowgrapgy 
and  Raman/  scattering  measurements  are  presented  here.  Metamorphic  behavior  of  the  jet  boundary  undergoing  disruption 
has  been  shown  to  occur  near  the  critical  point  of  the  injected  Ofluid.  At  subcritical  pressures,  the  formatioijjof  droplets  and 
ligaments  (fare  seen,  resembling  a  second  wind-induced  liquid-jet  breakup.  The  jet  is  incapable  of  reaching  the  classical  full 
atomization  regimeQdue  to  existence  of  the  critical  point  of  the  jet  fluid  Cwithin  the  operating  range.  At  this  point,  the  jet 
anatomy  changes  abruptly  to  imitate  turbulent  gas  jet  injection.  The  jet  initial  growth  rate  is  plotted  against  the  chamber-to- 
injectant  density  ratio,  along  with  available  data  on  other  liquid/gaseous  jets  and  mixing  layers,  producing  a  unique  and 
informative  graph.  For  supercritical  conditions,  our  measured  growth  rate  agrees  well  with  a  theoretical  equation  proposed 
by  Dimotakis  [1]  and  closely  follows  the  trend  of  Papamoschou  and  Roshko  [2]  for  incompressible  but  variable-density 
gaseous  turbulent  mixing  layers.  Fractal  analysis  of  the  jet  boundary  also  shows  a  similarity  to  gas  jet  behavior  with 
comparable  fractal  dimension.  This  is  the  first  time  quantitative  evidence  has  been  provided  to  support  qualitative 
visualizations  suggesting  that  supercritical  jets  appear  to  behave  like  conventional  gas  jets.  Based  on  a  physical  hypothesis 
implying  that  at  the  point  where  gas-like  behavior  appears,  the  value  of  the  jet  interface  bulge  separation/formation 
characteristic  time  should  be  near  that  of  gasification  time,  an  “intuitive/smart”  equation  is  proposed  dfor  the  growth  rate 
that  agrees  well  with  the  presented  experimental  data  for  a  wide  range  of  density  ratios.  Initial  Raman  scattering 
measurements,  being  proportional  to  density,  provided  similar  spread  rate  if  twice  the  FWHM  of  the  radial  density  profiles 
is  considered.  Radial  profiles  appear  to  show  similarity  behavior  but  hot  of  the  classical  type  observed  for  turbulent 
gaseous  jets  at  far  field. 


INTRODUCTION 


,7 


iiiiitvi/v/viiVM  a/  \/jr 

I  iP  t 

As  combustion  and  thrust  /  chamber  pressures  are 
raised  to  benefit  ©{/-performance  and/or  efficiency  in 
cryogenic  rocket  and  turbine  engines,  the  injected  liquid 
fuel  may  experience  an(ambienp  under  thermodynamic 
supercritical  condition.  As  examples,  the  Space  Shuttle 
main  engine  thrust  chamber  pressure  with  liquid  H2 
(Pc=1.27MPa,  Tc=33.25K)/liquid  02  (Pc=  5.043MPa, 
TC=154.58K)  is  about  22.3  MPa,  and  the  combustion 
chamber  pressure  for  Vulcain  (Ariane  5)  can  reach  up  to 
10  MPa.  The  uniqueness  of  the  situation  comes  from  the 
fact  that  for  Qpure  substances,  the  distinction  between 
liquid  and  gas  phases  disappears  at  and  above  the  critical 
pressure  and  0  density,  Qthermal  Oconductivity,£)mass 
diffusivity,  and  surface  tension  show  large  variations  near 
the  critical  point.  At  sufficiently  high  ambient  pressures, 
the  solubility  of  Ogases  into  the  liquid  phase  becomes 
important!)  and  0  multicomponent  phaseQequilibrium 
information,  or  “critical  mixing  temperature  or  pressure,” 
should  be  usecW(see  Bruno  and  Ely  [3^,  QThe  terms 
subcritical  and  ^Ipercritical,  and  the  reduced  pressure  (Pr) 
and  temperature  (Tr),  are  defined  here  with  respect  to  the 
critical  condition  of  the  pure  substance  used  in  the  Qjets, 
and  not  the  environment. 


Past  efforts  by  other  researchers  on  this  subject  can  be 
found  inflfer  example^  Newman  and  Brzustowski  [4]  and 


found  in^for  example^  Newman  and  Brzustowski  [4]  and 
Mayer  et  al[5],  [6],  (amp  Chen  and  Sui  [7]  at  high 
Reynolds  numbag^and  Woodward  and  Talley  [8]  at  low 
Reynolds  numbers.  In  summary,  they  reported  gas-jet  like 
visual  Qappearanc eO  under  supercritical  chamber 
temperature  and  pressure  with  no  evidence  of  droplet  even 
though  liquid  was  injected.  At  low  Reynolds  numbers  and 
for  an  N2-into-N2  cryogenic  jet,  it  was  shown  that  the 
addition  of  He  gas  to  the  supercritical  ambient  N2  fluid 
forced  the  jet  to  look  like  an  injection  into  a  subcritical 
ambient  condition.  0)Raman  scattering  was  used  by 
Anderson  et  al  [9]  ,  Decker  et  al.  [10],  Oschwald  and 
Schik.  [11],  and  Oschwald  et  al.  [12]  to  investigate  the 
structure  of  the  nitrogen  fluid  jet  under  supercritical 
condition.  On  these  experiments.  Decker  et  al  [10] 
observed  a  smooth  transition  from  the  high-density  core  of 
an  N2  into  N2  jet  toward  the  gaseous  outer  region;  that  is, 
there  was  no  distinct  phase  interface.  The  density  radial 
profile  was  never  a  top-hat  even  as  close  as  a  distance-to- 
diameter  ratio  of  2.5.  Based  on  limited  tests,  they 
concluded  that  the  thermodynamic  state  of  an  injected 
supercritical  fluid  is  of  prime  importance  in  the  jet  growth 
and  not  the  injection  velocity  or  momentum. 


The  objective  in  this  paper  is  to  concisely  present  our 
experiences  learned  so  far  on  injection  under  supercritical 
condition  through  experimental  and  modeling  efforts. 

DESCRIPTION  OF  THE  TEST  RIG 

J*' 

The  stainless  steel  chamber  in  Figo  1  can  withstand 
pressures  and  temperatures  of  up  to  137  atm  (2000  psi) 
and  473  K,  respectively,  and  has  two  facing  circular 
sapphire  windows  and  two  UV-grade  side-mounted  slot¬ 
shaped  quartz  windows  for  laser  in/out  of  the  chamber. 
Transport  of  the  liquid  through  the  lines  and  the  injector  is 
achieved  as  shown  in  Figure  1.  The  liquid  N2  from  a  dewar 
is  used  to  cool  and/or  liquefy  the  injectant  passing 
through  the  cryogenic  cooler  prior  to  injection.  In 
addition  to  N2,  a  branch  is  connected  for  possibly 
introducing  other  gases  or  mixtures  of  gases  as  a 
chamber  medium.  The  mass  flow  rate  of  the  injectant  is 
regulated  and  also  measured  via  a  mass  flowmeter,  and 
a  precision  micrometer  valve.  For  more  detailajefer  to 
Woodward  and  Talley  [8]  and  Chehroudi  et  al.  [13). 
The  injection  in  this  studyOis  through  a  sharp-edged  5fl© 
— (mn)  long  stainless  steel  tube  with  1.5<@pm 

diameter  and  a  25^iicron  (0.010!Z|,jgiier  hole  (length- 
to-diameter  ratio  of  200).  With  the  Reynolds  number 
ranging  from  25,000  to  70,000,  the  entrance  length  from 
50  to  100  is  neededj(see  Schlichting  [  14L The  length  is 
''therefore  long  enough  to  Oensure  mlly\developed 
turbulent  pipe  flow  at  the  exit.  The  rig  is  fully 
instrumented  with  thermocouples,  pressure  gages  and 
transducers,  and  mass  flowmeters  at  the  locations 
indicated  in  Fi&jL  Back-illumination  and  a  model  K2 
Infinity  long-distance  microscope  is  used  with  a  TM- 
745E  high  resolution  p68(H)  x  493(V)  pixels  in 
8.8(H)x6.6(V)  mm  actual  sensing  areaj  interlaced  CCD 
camera  by  PULNix  to  form  images  of  the  injected  jets. 


the  jet.  Major  structural/interface  changes  (Spccurred  at 
about  Prof  1.03.  No  drops  were  detected,  thread-  or  finger- 
like  entities  at  the  interface  were  not  broken  up  into 
droplets  but  dissolved  at  different  distances  from  the  jet 
central  dense-core  region.  This  created  a  mixing  layer  in 
which  phase  transition  and/or  large  local  density  non¬ 
uniformities  occurred.  A  further  increase  in  chamber 
pressure,  decreased  the  length  and  the  thickness  of  the 
central  dense-core  and  images  progressively  appeared 
similar  to  an  injection  of  a  gaseous  turbulent  jet  into  a 
gaseous  environment  (see  images  in  Fig.  2).  Aligned  with 
what  reported  before,  no  evidence  of  drop  formation  was 
seen  in  this  gas-like  jet  regime.  The  reason  for  this 
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For  more  detaijs^efer  to  Chehroudi  etal.  [13]. 
EXPERIMENTAL  AND  MODELING  EFFORTS 


In  this  section,  the  authoi^xperiences  learnei 


_ ned.at  the 

AFRI^are  described  and  summarized lio  this  dater  First, 
structural  changes  observed  using  shadowgraphy  are 
explained.  Second,  Qesults  of  fractal  analysis  is 
presented  and  briefly  discussed.  Third,  initial  attempt 
to  model  the  measured  growth  rate  of  the  jet  is  presented. 
Finally,  ew£  recent  effort  Qusing  laser  diagnostics  via 
Raman  scattering  0to  measure  density  profiles  and  growth 
rate  is  briefed.  This  will  provide  a  summary  of  mt. 
important  findings  and  achievements  regarding  injection 
under  supercritical  condition. 

VISUAL  CHANGES  OF  THE  JET  T 


)  Gaseous  N?  Ambient) 


Using  liquid  N2  injection  into  ^gaseous 
when  chamber  pressure  was  far  lower  than  the  critical 
pressure  of  the  jet  substance  (i.e^subcritical  chamber 
pressure),  the  jet  was  observed  (er  be  liquid-like  with 
surface  instabilities  that  grew  downstream  from  the 
injector.  Also,  very  fine  drops  were  detected  ejecting  from 


Figure  1.  Schematic  diagram  of  experimental  setup  for  sub*  to  supercritical  jet  injection 

metamorphic  behavior  is  to  be  sought  in  progressive 
reduction  of  surface  tension  and  heat  of  vaporization  until 
they  both  are  sufficiently  reduced  near  and  above  the 
critical  point  Similar  observations  were  made  for  the  02- 
into-N2  case  but  with  the  transition  to  gas-like  jet  atnear 
the  oxygen  critical  pressure  (Pr  =  0.85).  In  FigJ^ZT  as 
chamber  pressure  approaches  the  critical  condition,  surface 
tension  is  reduced  to  a  near-zero  value  and  for  the  range  of 
Reynolds  number  used  here  (25,000  to  75,000),  the 
Ohnesorge  number  sharply  swings  from  a  low  (estimated 
to  be  2.8xlO'3£)  to  a  very  large  (infinity  when  surface 
tension  is  zero)  values.  At  the  low  Ohnesorge  number  in 
this  work,  the  so-called  second-wind  induced  liquid  jet 
breakup  behavior  is  observed,  see  Fig^jT  Past  this  stage, 
and  before  reaching  the  full  atomization,  the  surface 
tension  is  sufficiently  and  rapidly  reduced  to  produce  a 
gas-like  jet  appearance  near  but  before  reaching  the  critical 
pressure  valu^(see  Fig.  ^  The  aforementioned  threads  or 
fingers  are  gasified  and  no  drops  are  detected.  Transition 
into  a  full  atomization  regime  is  therefore  inhibited  by  this 
mechanism. 
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QUANTITATIVE  MEASURE  OF  THE  JET 
GROWTH  RATE 

The  jet  growth  rates  or  spread  angles  were  measured 
from  the  images  and  only  to  within  a  ^gpn  distance  from 
the  injector  exitolane  (distance-to-diameter  ratio  of  up  to 
28).  Bated  oft  jhe  intact  core  measurements  of  the  liquid 
sprays  and  potential  core  of  the  gaseous  jet ^(Qhehroudi  et 
al  [13ty  showed  that  this  measurement  ^distance  was  well 
within  the  mixing  layer  region  and  that  comparison  with 
classical  mixing  layers  was  justifiei(see  Chehroudi  et  al 
[19]  and  Abramovich  [18]U  Detailed  discussion  of  the 
choices  of  the  data  set  in  Fi<£*3  is  given  elsewhere  in  our 
published  work^see  Chehroudi  etal.  [13fy, 

The  simplest  curves  in  Fi^3  are  the  prediction  of  the 
linear  jet  growth  for  turbulent  incompressible  submerged 
jets  using  the  mixing  length  concept  and  a  semi-empirical 
equation  by  Abramovich  [18]  incorporating  the  effects  of 
density  variations  using  a  characteristic  velocity.  The 
growth  rate  measurements  by  Brown  and  Roshko  [20]  in  a 
subsonic  two-dimensional  incompressible  turbulent  mixing 
layer  between  helium  and  nitrogen  are  shown  in  FigT3 
Also,  the  independently-derived  theoretical  equations 
proposed  by  Papamoschou  and  Roshko  [2]  and  Dimotakis 
[1]  for  incompressible  variable-density  mixing  layers  are 
shown. 


Because  of  the  liquid-like  and  gas-like  appearances  of 
the  jet,  the  growth  rate  for  liquid  sprays  from  single-hole 
nozzles  typical  of  the  ones  used  in  diesel  engines  are  also 
shown.  We  selected  the  steady  isothermal  spray-angle 
information  by  Reitz  and  Bracco  [15]  and  Hiroyasu  and 
Arai  [21]  at  two  different  nozzle  length-to-diameter  ratios 
with  their  corresponding  vertical  error-bands  indicating 
experimental  scatter.  In  addition,  a  theoretical  equation 
derived  by  Reitz  and  Bracco  [15]  is  shown.  As  a  cross 
check,  a  recent  curve-fitted  equation  to  experimental  data 
from  transient  sprays  proposed  by  Naber  and  Siebers  [22] 
is  also  shown.  Their  angle  measurement  zone  extends 
beyond  our  Qinitial  region  and  this  to  some  extent 
contributes  to  disagreement  seen  between  the  two  sets  of 
data  for  liquid  sprays  at  injector  length-to-diameter  density 
ratio  of  about  4,  see  Fi^x 

Figure  3  covers  a  density  ratio  of  nearly  four  orders  of 
magnitude,  from  liquid  sprays  to  supersonic  mixing  layers, 
a  unique  and  new  plot  on  its  own  right.  Disagreements 
seen  between  the  measured  data  can  to  some  extent  be 
attributed  to  differences  in  the  definition  of  the  mixing 
layer  thicknesses  and  their  measurement  method^(see 
Chehroudi  et  al.  [ 13j ,  The  constant  spread-angle  results  of 
incompressible  turbulent  jet  overpredicts  nearly  all  others 
in  Fi^.  Disagreement  between  the  turbulent  gas  jet  by 
Abramovich  [18]  and  incompressible  variable-density 
model  of  Papamoschou  and  Roshko  [2]  increases  as 
density  ratio  increases.  It  is  clear  that  for  a  range  of  density 
ratios  in  which  our  images  show  gas-jet  like  appearance, 
the  experimental  data  agrees  well  with  the  proposed 
theoretical  equation  by  Dimotakis  [1]  and  closely  follows 
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Figure  2.  (Top)  Back-illuminated  images  of  a  cryogenic 
N2  into  N2  round  jet  injection  at  chamber  conditions  of 
Pr  =  0.92,  1.22,  and  2.7  but  at  a  fixed  supercritical 
temperature  of  300  K  (PCriticai  =  3.39  MPa;  Tc  =  126.2  K). 
Re=  25,000  to  75,000,  Froud  number:  40,000  to  110,000 
(momentum-dominated  jet),  injection  velocity:  10-15  m/s, 
and  injectant  temperature:  99  to  110  K.  (Bottom) 
Software-magnified  images  of  the  edge  of  the  jet  in 
transition  to  the  gas-jet  like  appearance  is  seen  starting  at 
just  below  the  critical  pressure  of  the  injectant.  A  gradual 
transition  from  a  classical  liquid-like  appearance  in  which 
ligaments  and  drops  are  formed  at  the  boundary  in  the 
liquid  atomization  regime  to  a  comb-like  structure  near 
the  critical  point,  and  finally  to  where  a  submerged 
turbulent  gas  jet  appearance  emerges  can  be  observed. 


the  trend  of  Papamoschou  and  Roshko  [2],  This  was  the 
first  quantitative  confirmation  that  «&n  injected  jet  under 
supercritical  pressure  and  temperature  ambient  conditions 
(based  on  the  pure  injectant)*  grows  similar  to  a  gas  jet. 
Even  though  the  jet  investigated  here  appears  to  go 
through  initial  phases  of  the  liquid  atomization  process, 
there  is  a  marked  disagreement  at  above  the  critical 
condition  in  both  the  magnitude  and  slope  between  the 
liquid  sprays  (at  a  comparable  length-to-diameter  ratio  of 
85)  and  our  data  in  Fij£3.  It  appears  that  although  the  jet 
studied  here  shows  second  wind-induced  breakup  features 
similar  to  liquid  jets,  it  fails  to  reach  full  atomization  state 
as  chamber  pressure  (really  density)  is  raised.  This  is 
because  the  thermodynamic  stke  approaches  the  critical 
point  and  consequently  both  surface  tension  and  heat  of 
vaporization  are  reduced  to  sufficiently  low  values. 
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Jet  or  Spray  or  Mixing  Layer  Divergence  Angie  vs  Chamber-to-injectant 
_ Density  Ratio _ 
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Figure  3.  Shows  spreading  or  growth  rate  as  a  tangent  of  the  visual  spreading  angle  versus  the  chamber-to-injectant  density  ratio.  (*) 
refers  to  data  taken  at  AFRL.  —  'ut 

FRACTAL  ANALYSIS  OF  THE  JET 

This  appeared  to  be  a  suitable  analysis  method  to 
investigate  the  morphology  of  the  interfacial  phenomena 
and  in  recent  years  a  number  of  applications  of  fractal 
analysis  have  been  demonstrated  in  different  disciplines. 

For  a  concise  description  and  methodology  refer  to 
Chehroudi  et  al.  [23].  For  example,  Sreenivasan  and 
Meneveau  [24]  computed  the  so-called  fractal  dimension 
of  the  turbulent/non-turbulent  boundary  of  an 
incompressible  axisymmetric  gaseous  jet  and  reported  a 
value  of  1.33  using  two-dimensional  slicing  by  a  laser 
sheet.  Sreenivasan  [24]  also  reported  values  of  1.35^1.34, 
and  1.38  for  a  round  gaseous  jet,  a  plane  gaseous  mixing 
layer,  and  a  boundary  layer  flow,  respectively.  Figure  4 
shows  results  of  plotting  the  fractal  dimension  applied  to 
the  N2-into-N2  jets  as  a  function  of  relative  chamber 
pressure.  In  this  figure  the  fractal  dimension  by  other 
researchers  in  the  gaseous  jets,  mixing  layer,  and  boundary 
layers  are  also  shown.  The  key  observation  here  is  that  the 
average  value  measured  for  our  jet  reaches  very  near  those 
of  gaseous  jets  and  mixing  layers  when  ambient  pressure 
goes  above  the  critical  point  of  the  injectant.  This  is 
considered  as  an  additional  quantitative  confirmation  of 
the  gas-like  jet  behavior  indicated  earlier  through  our 
growth  rate  measurements  in  this  same  region.  As 
pressure  (density)  is  decreased  below  Pr  of  about  0.6,  the 
fractal  dimension  is  rapidly  reduced  towards  the  Euclidean 
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Figure  4.  Box-counting  and  Minkowski  (Euclidean 
Distance  Mapping  (EDM)  algorithm)  fractal  dimensions  of 
the  visual  boundary  of  the  jet  as  a  function  of  the  reduced 
chamber  pressure  for  N2-into-N2  injection. 


» 


value  of  1  for  a  smooth  circular  cylinder  with  no  surface 
irregularities.  Considering  the  geometric  interpretation  of 
the  fractal  dimension,  measured  values  are  quite 
consistent  with  the  visual  observations  of  jet  interface 
changes  in  our  tests.  For  more  details  see  Chehroudi  et  al. 
[23).  ' 

MODELING  OF  THE  GROWTH  RATE 

As  no  analytical  treatment  exists  starting  from  basic 
principles,  the  idea  was  to  intelligently  search  for  a  form 
of  an  equation  based  on  an  intuitive  and  physically 
reasonable  model  tiiat  can  predict  the  jet  growth  rate  data 
presented  in  FigjTL  In  essence,  the  proposed  approach  is 
to  employ  two  appropriate  characteristic  times  of  the 
problem,  combined  with  a  hypothesis  and  the  growth  rate 
equations  for  isothermal  0  liquid  sprays  and  Qgaseous 
mixing  layers  to  propose  an  “intuitive/smart”  equation 
that  is  capable  to  closely  represent  the  growth  rate 
information  measured  for  the  sub-  and  supercritical 
conditions  presented  before. 


and  drops  because  they  are  gasified  as  fast  as  they  desire  to 
be  detached.  In  a  sense,  this  is  defined  as  the  onset  of  the 
gas-jet  like  behavior. 

It  was  indicted  earlier  that  the  0jets  investigated  here 
showed  both  liquid-like  and  gas-like  behaviors  and  that  the 
growth  rate  of  the  latter  agrees  well  with  those  of 
incompressible  variable-density  gaseous  mixing  layers. 
Therefore,  forms  of  the  growth  rate  equations  for  the 
liquid  sprays  as  well  as  incompressible  variable-density 
gaseous  jets  are  considered.  For  liquid  sprays,  using  linear 
instability  theoretical  analysis,  Reitz  and  Bracco  [15] 
proposed  the  following  equation  for  the  isothermal  liquid 
spray  growth  rate:  0ss  0.27  [0  +  (pg/pi)0'5],  where  0  is 
spray  Q  full-cone  angle.  £^?Also,  Qfrom 
Brown/Papamoschou/Rashko’s  Otheory  O for  O  an 
incompressible  variable-density  gaseous  mixing  layer  one 
can  derive  an  equation,  0B/p/rS  0.212  [1  +  (pg/pi)0  5],  where 
0  is  the  visual  jet  full-cone  angle.  Finally,  based  on  a 
different  theoretical  treatment  by  Dimotakis  for 
incompressible  variable-density  gaseous  mixing  layers,  an 


Consider  injection 
into  a  subcritical  pressure 
environment  similar  to  die 
ones  shown  in  Fig*2. 

Here,  an  appropriate 
characteristic  time  is  the 
“bulge” 

formation/separation  time 
(tb)  pertaining  to  the 
interface  of  turbulent 
liquid  jet  where  isolated 
ligaments  and  drops  are 
produced  through 

primary  atomization. 

Tseng  et  al.  [26] 
suggested  that 

xb=(p]L3/o)1/2  for  the 
primary  breakup  of 
turbulent  liquid  jets. 

Here,  pj,  L,  and  a  are 
liquid  density, 

characteristic  dimension 
of  turbulent  eddy,  and 
surface  tension, 

respectively.  The  second 
characteristic  time  (at 

subcritical  pressures)  is  Figure  5.  Comparison  of  the  proposed  growth  rate  model  with  the  experimental  data  for 
the  gasification  time  (Xg)  the  sub-  to  supercritical  jet.  (*)  refers  to  data  taken  at  AFRL. 
calculated  through  the  so- 


called  D-square  law  for 

drops  and  is  equal  to  D2/K  where  D  and  K  are  drop 
diameter  and  vaporization  constant,  respectively.  In 
addition,  we  also  propose  a  hypothesis  that  if  the 
aforementioned  characteristic  times  (calculated  for 
appropriate  length  scales)  are  nearly  equal  in  magnitude, 
then  the  interface  bulges  are  not  able  to  separate  as 
unattached  entities  from  the  jet  interface  to  form  ligaments 


equation  of  the  form  0d=  0.265  [0.59  +  (Pg/pi)0'5]  can  be 
derived.  The  similarity  of  these  equations  is  remarkable 
despite  their  different  theoretical  treatments.  The  following 
“intuitive/smart”  form  of  an  equation  describing  the 
growth  rate  for  both  sub-  and  supercritical  jets  is  therefore 
proposed:  0^=  0.27  l('Zb/(xb+  Tg))  +  (pg/p^^see 
Chehroudi  et  al.  [23]  for  more  details).  In  a  sens&j-a 
physical  meaning  is  attributed  to  the  fust  term  in  me 


bracket  and  it  is  linked  to  the  physical  picture 
hypothesized  earlier.  Note  that,  in  the  limit,  this  equation 
collapses  to  the  isothermal  liquid  spray  case  when  Tg  >>  Tb 
and  Tg  oo. 

— —  Figure  5  shows  comparison  of  the  proposed  model 
equation  with  the  characteristic  time  ratio  calculated  for 
HP  ft/  the  N2-into-N2  Case.  Good  agreement  is  observed.  For 
cases  other  than  the  N2-into-N2,  such  as  N2-into-He  or  N2- 
into-Ar,  due  to  lack  of  reliable  information  on 
thermodynamic  properties  and  in  particular  the  surface 
tension.  To  model  these  cases,  it  was  proposed  that  the 
,  characteristic  time  ratio  Tb/(Tb+  Tg)  is  a  dominant  function 

™  of  the  density  ratio  only,  i.e.  Tb/(Tb+  Tg)  =F(Pg/pi).  As 

Brown  and  Roshko  [20]  indicated,  this  is  a  reasonable 
behavior  to  consider  because  at  low  Mach  number  squared 
values  there  is  no  distinction  between  mixing  layers  made 
of  two  streams  with  differences  in  molecular  weights, 
differences  in  temperature,  and  high  speed  compressibility 
effecty(see  also  Chehroudi  et  al.  [13p,Il  was  shown  that 
the  same  function  F(Pg/pi)  calculated  through 
measurements  for  the  N2-into-N2  case  could  be  made  to 
work  for  other  cases  provided  a  case-dependant 
transformation  is  made  to  the  density  ratio  at  which  F  is 
evaluated.  .  The  quality  of  the  agreement  with 
experimental  data  is  demonstrated  in  Fig£o  and  considered 
to  be  good. 

RAMAN  SCATTERING  MEASUREMENTS 


Figure  6.  Radial  plots  of  (I-  Ich),  the  Raman 
intensity  (I)  from  the  jet  as  measured  with  reference 
to  that  of  the  chamber  gas  (Ich)  for  the  N2  into  N2 
case  under  sub-  and  supercritical  conditions  and  at 
two  different  axial  positions. 


When  a  transparent  medium  is  irradiated,  some 
fraction  of  tjje^beam  is  scattered  in  all  directions.  Ramanr 
ju  1928,*oiscovered  that  the  wavelength  of  a  very  small 
fraction  of  the  scattered  radiation  by  certain  molecules 
shifted  in  an  amount  dependent  upon  the  chemical 
structure  of  the  molecules  responsible  for  the  scattering.  It 
can  be  shown  that  the  intensity  of  the  signal  is  proportional 
to  the  density  of  the  observed  molecule.  In  our  case,  an 
Nd:YAG  pulsed  laser  beam  at  532  nm  was  used  with  a 
passive  pulse  stretcher  to  change  the  pulse  duration  from 
10  ns  to  aboutO  36  ns.  This  was  done  to  minimize  the 
occurrence  of  the  stimulated  Raman  effects  particularly  at 
sub-  and  near-  critical  conditior^see  Anderson  et  al  [9] 
The  transmitting  optics  created  a  Ocontrollable  sheet  of 
laser  with  precise  traversing  capabilities  in  two  directions. 
The  scattered  light  at  607  nm  is  detected  by  a  256?g56 
CCD  camera  equipped  with  suitably-chosen  number  of*^ 
filters  to  maximize  the  signal-to-noise  ratio.  For  more  i_ 
details  on  experimental  setup  and  data  acquisition  refer  to 
Chehroudi  et  al.[27). 

The  integrity  of  the  entire  Raman  measurement 
system  was  tested  by  measurement  of  the  density  of  a  pure 
nitrogen  gas  in  the  high-pressure  chamber  for  several 
pressures.  Measured  to  calculated  (based  on  ideal  gas 
law)  density  ratios  agreed  to  within  2%.  Figure  6  shows 
intensity  defect,  being  proportional  to  the  density  defect, 
radial  plots  at  sub-  and  supercritical  conditions.  It  is  clear 
that  as  chamber  pressure  is  increased  the  intensity  level 
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Figure  7.  Normalized  plots  to  inspect  similarity  behavior 
in  comparison  to  gaseous  jet  results.  Laser  enters  from 
right  to  left  N2  into  N2  at  supercritical  chamber  condition. 

drops  down  in  the  central  region  due  to  larger  spreading  of 
the  jet.  The  data  for  the  subcritical  case  should  be 
considered  with  reservation  particularly  at  larger  radial 
positions  due  to  reflection  of  the  laser  from  the  jet’s  well- 
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Figure  8.  Plot  of  the  spread  rate  as  measured  by  Raman  intensity  Dradial  plots  atOtwice  the  FWHMQshows 
agreement  with  our  earlier  measurements  using  shadowgraphy  technique. 


defined  boundary  causing  broadening  of  the  profiles.  This 
issue  is  addressed  in  more  details  in  Chehroudi  et  al  [27]. 
If  the  jet  exhibits  gas-like  behavior  one  might  attempt  to 
compare  the  results  with  those  of  gaseous  jets  to  see  extent 
to  which  they  simulate  each  other.  A  normalized  plot 
format  is  used  for  this  purpose  in  Figure  7.  Here,  results  at 
supercritical  condition  and  for  different  axial  positions 
are  shown  along  with  a  similarity-type  model  equation 
generally  accepted  for  turbulent  jets  andat  larger  axial 
distances  than  those  shown  in  thts-pj'ure^Even  though? 
some  degree  of  similarity  is  exhibited,  they  are  not  of  the 
type  observed  in  the  far  field  region  for  turbulent  gaseous 
jets.  Actual  density  measurements  in  gaseous  jets  made  by 
So  et  al.  [28]  in  comparable  nondimensionalized  axial 
distances  indicate  reasonable  agreement  within  the 
experimental  accuracy^see  Fig.  ly  Results  by  Oschwald 
et  al.  [12]  are  also  shown  but  no  conclusion  is  drawn  at 
this  time  due  to  differences  in  initial  conditions  between 
the  two  sets  of  data  known  to  affect  jets  and  sprays 
spreading.  The  full  width  of  the  radial  profiles  measured 
at  the  half  maximum  (FWHM)  is  extracted  fronrlhe 
Raman  measured  data.  They  are  believed  to  be  unaffected 
by  the  broadening  observed  at  the  edge  of  these  profile. 
When  spread  rate  is  calculated  based  on  twice  the  FWHM 
values  of  Raman  radial  profiles  and  plotted  along  with 
other  measurements  in  Fig^B,  it  shows  good  agreement 
with  our  earlier  results  using  shadowgraphy  method.  This 
can  be  considered  as  a  correspondence  calibration  of  the 
two  different  techniques  of  estimating  the  width  of  the  jet 
investigated  here. 

SUMMARY  AND  CONCLUSIONS 

Metamorphic  behavior  and  growth  rates  of  jets 
injected  into  an  environment  at  a  given  supercritical 


temperature  but  varying  pressure  are  analyzed.  At. 
subcritical  chamber  pressures,  the  injected  fluid 
appears  to  undergo  the  classical  liquid  jet  breakup 
stages  up  to  a  second  wind-induced  breakup  regime 
with  discemable  ligaments  and  droplets.  Penetration 
into  the  full  atomization  regime  is  inhibited  near  the 
critical  pressure  of  the  injectant  because  of  the 
combined  effects  of  lowered  surface  tension  and  heat 
of  vaporization.  At  this  point  the  visual  structure  of 
the  changes  to  a  gas-jet  like  appearance  and  remains 
up  to  the  highest  pressure  tested.  Also,  a  unique  and 
comprehensive  plot  on  growth  rate  is  formed  using  the 
most  relevant  works  of  other  researchers  covering  an 
ambient-to-injectant  density  ratio  range  of  four  orders  of 
magnitude.  Our  measured  jet  growth  rates  via 
shadowgraphy  follow  the  theoretical  equations  proposed 
by  Papamoschou  and  Roshko  [2]  and  Dimotakis  [1],  both 
derived  for  incompressible  variable-density  gaseous 
turbulent  mixing  layer.  This  agreement  starts  at  a  pressure 
near  but  below  the  thermodynamic  critical  pressure  of  the 
injectant,  quantitatively  supporting  the  observed  gas-jet 
like  visual  appearance  of  the  supercritical  jets  for  the  first 
time.  The  geometry  of  the  jet  interface  has  also  been 
examined  for  the  first  time  by  fractal  analysis.  The  results 
clearly  indicate  a  transition  from  a  Euclidean  to  a  fractal 
interface,  with  a  fractal  dimension  close  to  values 
measured  for  gaseous  turbulent  jets.  This  provides 
additional  quantitative  evidence  for  the  hypothesis  that  the 
jet  evolves  into  a  gas-like  behavior.  An  “intuitive/smart” 
equation  is  proposed  that  agrees  well  with  the 
experimental  growth  rate  data,  based  on  a  proposed 
physical  mechanism  and  characteristic  gasification  times 
and  interfacial  bulge  formation/separation  times.  Initial 
Raman  scattering  measurements  from  the  jets  confirms  the 


growth  rate  measured  via  shadowgraphy  provided  ftwice 
the  FWHM  of  the  ^density  radial  profiles  are  used.  The 
density  defect  plots  show  tendency  to  be  similar  but  not  of 
the  type  observed  for  turbulent  gaseous  jets  at  far  field. 
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